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The non-Markovian dynamics of quantum entanglement is studied by the Shabani-Lidar master 
equation when one of entangled quantum systems is coupled to a local reservoir with memory effects. 
The completely positive reduced dynamical map can be constructed in the Kraus representation. 
Quantum entanglement decays more slowly in the non-Markovian environment. The decoherence 
time for quantum entanglement can be markedly increased by the change of the memory kernel. 
It is found out that the entanglement sudden death between quantum systems and entanglement 
sudden birth between the system and reservoir occur at different instants. 
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1. Introduction 

The quantum entanglement in composite systems is a key issue in quantum information 
theory [l|, The resources of quantum entanglement are contained in nonclassical states 
which are useful for quantum information processing, like quantum teleportation, cryptog- 



raphy, dense coding and quantum computation 



It is known that entanglement can 



be efficiently quantified by some useful measurement including the concurrence [6j, the neg- 

n n n 

ativity [7J and relative entropy [8|, |9[. Because of interactions between open systems and 
surrounding environments, the decoherence of quantum entanglement is usually inevitable 
101 ] . Therefore, dynamical properties of quantum entanglement have received much atten- 



tion 



llH24| . In some works 



25 



26] , the dynamics of the two-qubit entanglement have been 



investigated when the noninteracting qubits are respectively coupled to two independent 
reservoirs. It was found out that the sudden death and revival of the entanglement can 
happen at some instants. The Markovian approximation is used for the study of dynamics 
of open systems weakly coupled to reservoirs without memory effects. The corresponding 
map is always completely positive and trace-preserving. However, in realistic environments 
with memory effects, time evolutions of quantum systems often obey the non-Markovian 



dynamics 



10| . The memory effects can be described by some phenomenological master 



equations with a reasonable memory kernel 27H3l|. It is an obstacle to deal with physical 
decoherence of open systems analytically and numerically in non-Markovian environments. 
Among these master equations, the Shabani-Lidar one interpolates between the exact dy- 
namics and Markovian one 29|. The desirable property of Shabani-Lidar master equation 



is the preservation of the complete positivity of the reduced dynamics. The requirement of 
complete positivity of dynamics ensures that the physical states of reduced systems can be 
applied to practical tasks of quantum information processing. From the point of view of a 
definite physical process, it is of great value to obtain a completely positive dynamical map 
in a non-Markovian bath and investigate the dynamics of quantum entanglement. 

In this paper, by the Shabani-Lidar master equation, we construct the completely-positive 
dynamical map in order to analyze the decoherence of a two-level quantum system coupled to 
a vacuum reservoir with general exponential memory. In Sec. 2, the reduced dynamical map 
is written in the Kraus representation. The dynamics of quantum entanglement between 
the open qubit and an ancillary system isolated from the local reservoir is obtained in Sec. 
3. The memory effects on distribution of quantum entanglement between the systems and 
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local bath are also taken into account. A simple discussion concludes the paper. 



2. The completely positive reduced dynamical map 

The decoherence of quantum system is always unavoidable because of the interactions 
with its local reservoir. To set up the reduced dynamical map, we investigate the total 
Hamiltonian of the effective two- level open system and the local environment, 

H = H S + H E + H SE , (1) 

where 

H s = u a + a_ (2) 
describes the intrinsic Hamiltonian of the open two-level system, 

H E = J2 UJ j a ] a o ( 3 ) 

i 

is the Hamiltonian of the reservoir with a large amount of independent harmonic oscillators 
and 

H S e = ^2 9j(v+aj + a -a)) (4) 

represents the interaction between the open qubit S and the environment E. The parameter 
u>o describes the transition energy from the ground state \g) to the excited one |e). The 
operators a± = (a x ± icr y )/2 are the raising and lowing operators of the open qubit. Here 
o)j , cij are the creation and annihilation operators of the reservoir. In the environment with 
memory effects, the time evolution of an effective two-level system can be approximately 
obtained by the Shabani-Lidar master equation [29]. The master equation can give new 
insight into the dynamics of non-Markovian quantum systems. The reduced states for the 
open system ps{t) can be written as 



dps{t) 
dt 



C [ f(t')exp(a')p s (t-t')dt'. (5) 
J o 



Here the rotating-wave approximation is adopted. In the following analysis, the initial total 
state ps(0) <g> |0)b(0| is assumed. \0)e denotes the vacuum state of the local reservoir. The 
fucntion f(f) is the phenomenological kernel characterizing the memory effects from the 
environment. The memory kernel function can be experimentally determined by quantum 



state tomography 29j. In order to understand the physical meaning of the Shabani-Lidar 
memory kernel, we stress that this master equation describes a situation in which a mea- 
surement of the environment at time t' is followed by a Markovian evolution, described in 



terms of continuous measurements of the environment at times t > t' |29|, |3jJ. The time 
t' characterizes the bath memory effects. The memory kernel is a function which assigns 
weights to different measurements. The determination of the memory kernel is related to 



the reservoir spectral density. The local Liouvillian operator C can be expressed by 



(6) 



The parameter 70 = 7o(<?j, Wj, coo) describes the decaying rate in the Markovian model. 
For the dynamical map $(£) of the reduced system S, the density matrix at any time t 
can also be expressed by psif) = ®(t)ps(0). According to [29], the definite expression 
of the dynamical map is given by $(t)ps(0) = Y^=o OWT r [£jPs(0)]-R;. Here Lj and Rj 

n 

are the j-th left and right damping basis of the Liouvillian operator [321 ] . In this case, 
£j(t) = Lap _1 [- — ^ Tpjj} — )^~y] w h ere {Aj = 1, —70, —7o/2, —70/2} are the four eigenvalues 
of the Liouvillian operator and the Laplace transformation F(s) = Lap[/ (£')]. 

To investigate the conditions of complete positivity of the dynamics, the Choi matrix V 
need be constructed by the elements V(i,j) = $>(t)\i)(j\ (33]. In the context, the general 
exponential memory kernel function is considered as f(t) = v4exp(— 7^) where A is the 
amplitude of the kernel and 7 represents the memory decaying rate. This kind of memory 
kernel can be physically modeled by 30|. The definite expression of the matrix V can be 
obtained by 

/ £(7o,0 £(7o/2,t)\ 

l-£(7o,t) 


V^(7o/2,t) 1/ 

where the elements £,(x,t) can be written as £(x,t) = exp[— ( - 1+ ^' |a:t ](cosh xut + sinhxtit). 
The parameter is u = | a/(1 + g) 2 — 4a|/2 where a = A/x and g = j/x. If a > (1 + (7) 2 /4, 
the above function £(x, t) is obtained by substituting sinh[-] and cosh[-] with sin[-] and cos[-]. 

According to the results of [29], the complete positivity of the map is equivalent to the 
positivity of the matrix V where all of the eigenvalues are non-negative. The condition of 



V 



(7) 
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complete positivity of the map can be analytically given by the inequality of 



£ 2 (7o/2,t)<£( 7 o,t)<l. (8) 

The impacts of the memory kernel on the complete positivity of the dynamics are numerically 
calculated in Fig. 1. When A < 7 or a < g, the values of D — £(70, t) — £ 2 (7o/2,£) are 
always non-negative in Fig. 1(a). Because the values of £(x, t) are monotonically decreasing 
in time, the inequality of £(x, 00) = < £(x, t) < £(x, 0) = 1 is always satisfied. This means 
that the corresponding map at any time is completely positive in the conditions of a < g. 
However, if the amplitude of the memory kernel satisfies a > (1 + g) 2 /4, the difference D 
oscillates between the positive values and negative ones in Fig. 1(b). For the time interval 
7o£ G (it — arctan 2n — arctan j^), the values of £(70, t) are negative which violates the 
positivity of the map. The region of D < represents the non-Markovian dynamical map 



which contradicts the complete positivity. According to the results of 311], during the time 
evolution which violates the positivity condition, the density matrix looses the probablistic 
interpretation. The reduced states in this regions are not physical. Therefore, the small 
amplitude of the kernel function is useful for setting up the completely-positive dynamical 
map. In the following analysis, we consider the physical evolution of the reduced states in 
the region of a < g. 

In the condition of the complete positivity, the reduced dynamical map can be expressed 
in the Kraus representation, 

Ps(t) = J2M k p s (0)Ml (9) 

k 

The Kraus operators {Mk} can be deduced by the completely-positive master equation with 
memory kernel. By the method of [29j , the ith column of the Kraus operator is the ith 
segment of the kth eigenvector of the completely-positive matrix V. The Kraus operators 
are written as 

M = y/m(a + \e) s (e\ + b+\g) s (g\) 
Mi = V±-^o,t)\g) s (e\ 

M 2 = y /jn{a.\e) s {e\ + b-\g) s (g\), (10) 

where the parameters /i ± = {1 +£(70, t) ± a/[1 - £(70, t)} 2 + 4£ 2 (7 /2, t)}/2, a±/b± = (/x ± - 
l)/£(7o/2,t) and b± = £(70/2, t)/ v / £ 2 (7 /2, t) + (/x± - l) 2 . For the Markovian reservoir 



without memory, the Kraus operators can be simplified into M = He)s(e| + \g)s(g\ an d 
Mi = y/\ — u 2 \g) s (e\ where v{t) = exp(— j t/2). 



3. Dynamics of quantum entanglement 

As is known, decoherence of quantum entanglement actually happens when local mea- 
surements are applied in quantum information processing 0, loj]. A reasonable model is 
considered as one open qubit S entangled with an ancillary qubit A which is isolated from 
a local vacuum reservoir E jjyj]. The decoheren ce p rocess can also be obtained by the map 



in terms of the total system-environment state 35 
total system can be obtained by 



36[. The evolution of the states for the 



p(t) = U SE (t)\p SA (0) <g> \o) E {o\]ul E (t), 



(11) 



where psa(O) is the initial entangled state. The unitary operator UsEif) can be expressed 
by the Kraus operators, 



Us E (t)\m)s®\0) E =J2Mk\m)s®\k)i 



(12) 



To clearly describe the dynamics of quantum entanglement between the system and local 
reservoir, we consider the negativity N |7| as one general measure for quantum entanglement. 
Based on the separability principle, the negativity N is introduced by, 



N(p) = max{0,-2^ ei } 



(13) 



where 6j is the ith negative eigenvalue of the partial transpose of the mixed state. For the 
separability of unentangled states, the partial transpose matrix has nonnegative eigenvalues 
if the corresponding state is unentangled. 

When the maximally initial state Psa(0) = |(l ee ) + \99))(( ee \ + (<7<?|) is chosen, the 
reduced state psEif) = Tr^[p(t)] at any time is written in the Hilbert space spanned by 
{|ee), \eg), \ge), \gg)}, 

I k(7o,t) |£(7o/2,*)\ 



PSA(t) 
















l-g(70,0 



V^(7o/2,t)0 







l 

2 



(14) 



The decaying of quantum entanglement N S a is shown in Fig. 2(a). It is seen that the 
values of N S a are always decreased to zero. The effects of memory kernel function on 
the entanglement are considered. When the relative amplitude a of the memory kernel 
f(t) is decreased, the vanishing of the entanglement N$a occurs more slowly. To clearly 
describe the decaying process, we can define the decoherence time r where the value of the 
entanglement is declined to 1/e. Fig. 2(b) clearly shows that the decoherence time scale 
is largely decreased with the increasing of the relative amplitude. We demonstrate that 
the decaying of quantum entanglement is much slower in the reservoir with the memory 
effects. This point is beneficial for the implementation of quantum computation in the 
non-Markovian environments. 

It is also interesting to study the distribution of quantum entanglement between the 
systems and local reservoir. After the tracing of the ancillary system A, the reduced states 
between the open system S and the local environment pse can be expressed in the Hilbert 
space spanned by {|e0), |el), \e2), \g0), \gl), \g2)}, 
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(15) 



where the elements are obtained by u = /i+a+, v = w = ^//i+^-a+a., z = 

yJp+P-b+b^ and q± = yfjl±a± ^/l — | (70, t) . The negativity iVs£ is also calculated and 
plotted as functions of the relative amplitude a and the time scale 7ot in Fig. 3(a). It is 
shown that the sudden birth of entanglement occurs at some time in the decoherence process 
and then the values are rapidly decreased to zero after enough time. With the decreasing 
of a, the time interval for the increasing of Nse is enlarged. Meanwhile, we also investigate 
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the dynamics of the reduced states pae after the tracing of the open system S, 
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(16) 



where some elements are given by p± = ,/p±a±. By the numerical calculation, the values of 
the entanglement Nae are always increased to the maximum after enough time in Fig. 3(b). 
This means that the initial entanglement can be completely transferred to that between the 
ancillary system and local bath during the decoherence process. It is also found out that 
the entanglement between the ancillary system and independent bath is produced when the 
open system is coupled to the local reservoir. 



4. Discussion 

We investigate the decoherence model of the open two-level system coupled to a local 
vacuum reservoir using the Shabani-Lidar master equation with memory kernel. For the 
general exponential memory, the small amplitude of the kernel function and the large decay- 
ing rate are helpful to construct the completely-positive reduced map. To clearly describe 
the dynamics of quantum entanglement, we construct the efficient map in the Kraus rep- 
resentation. By means of the negativity, the entanglement distribution between the system 
and local environment is also studied. It is found out that the decreasing of quantum entan- 
glement can be delayed pronouncedly in the non-Markovian reservoirs with memory effects. 
From the point of view of practical quantum information processing, the phenomenon is 
helpful for the understanding of the non-Markovian quantum computation. During the de- 
coherence process, the sudden birth of entanglement between the open system and local 
bath happens at some time. After enough time, the initial entanglement will be transferred 
to that between the ancillary system and independent environment. 
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Figure caption 

Figure 1 

(a). The difference D characterizing the complete positivity of the map is plotted as a 
function of a and 7ot when the parameters are chosen by a < g = 0.5; (b). The difference 
D is plotted when the parameters a >> g — 0.1. 

Figure 2 

(a). The dynamics of the entanglement between the open system and ancillary one N S a 
is plotted when g = 0.5. The solid line denotes the case of a — 0.1, the dashed one 
represents that of a — 0.3 and the dotted line describes the case of Markovian reservoirs; 
(b). The decoherence time scale 70T is numerically calculated with the change of the relative 
amplitude a of the kernel function. 

Figure 3 

(a). The entanglement N S e between the open system and local reservoir is plotted as 
functions of 7ot and a when g = 0.5; (b) The entanglement Nae between the ancillary 
system and reservoir is also plotted. 
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